Soliton solutions in an effective action for SU(2) Yang-Mills theory: including effects 

of higher-derivative term 
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The Skyrme-Faddeev-Niemi (SFN) model which is an 0(3) a model in three dimensional space 
upto fourth-order in the first derivative is regarded as a low-energy effective theory of SU(2) Yang- 
Mills theory. One can show from the Wilsonian renormalization group argument that the effective 
action of Yang-Mills theory recovers the SFN in the infrared region. However, the thoery contains 
an additional fourth-order term which destabilizes the soliton solution. In this paper, we derive 
the second derivative term perturbatively and show that the SFN model with the second derivative 
term possesses soliton solutions. 
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I. INTRODUCTION 

The Skyrme-Faddeev-Niemi (SFN) model which is an 
0(3) cr model in three dimensional space upto fourth- 
order in the first derivative has topological soliton so- 
lutions with torus or knot-like structure. The model 
was initiated in 70 's and its interest has been ex- 
tensively growing. The numerical simulations were per- 
formed in Refs. II II mi 01 > the integrability was shown 
in Ref. Ml, and the application to the condensed matter 
physics [3 and the Weinberg-Salam model ^ were also 
considered. The recent research especially focuses on the 
consistency between the SFN and fundamental theories 
such as QCD 113, IH El El. In those references, it 
is claimed that the SFN action should be deduced from 
the SU(2) Yang-Mills (YM) action at low energies. One 
can also show from the Wilsonian renormalization group 
argument that the effective action of Yang-Mills theory 
recovers the SFN in the infrared region However, 
the derivative expansion for slowly varying fields n upto 
quartic order brings an additional fourth-order term in 
the SFN model to destabilize the soliton solution. 

Similar situations can be seen also in various topolog- 
ical soliton models. In the Skyrme model, the chirally 
invariant lagrangian with quarks produces fourth order 
terms after the derivative expansion and they destabilize 
the soliton solution IB', To recover the stability of 
the skyrmion, the author of Ref. 17] introduced a large 
number of higher order terms in the first derivative whose 
coefficients were determined from the coefficients of the 
Skyrme model by using the recursion relations. Alterna- 
tively, in Ref. [lj| Gies pointed out the possibility that 
the second derivative order term can work as a stabilizer 
for the soliton. 

In this paper, we examine the Gies's supposition by 
numerical analysis. In section^ we give an introduction 



to the Skyrme-Faddeev-Niemi model with its topological 
property. In section ITTTl we show how to derive the SFN 
model action from the SU(2) Yang-Mills theory. In sec- 
tion llVI soliton solutions of this truncated YM action are 
studied. In order to find stable soliton solutions, we intro- 
duce a second derivative term which can be derived in a 
perturbative manner. The naive extremization scheme, 
however, produce the fourth order differential equation 
and the model has no stable soliton solution. Failure of 
finding the soliton is caused by the basic feature of the 
second derivative field theory. In section the higher 
derivative theory and Ostrogradski's formulation are re- 
viewed. We show the absence of bound state in the sec- 
ond derivative theory using an example in quantum me- 
chanics and introduce the perturbative treatment for the 
second derivative theory. In section IVll we present our 
numerical results. In section lVlIl are concluding remarks. 



II. SKYRME-FADDEEV-NIEMI MODEL 

The Faddeev-Niemi conjecture for the low-energy 
model of SU(2) Yang-Mills theory is expressed by fol- 
lowing effective action: 



^^{d^nf + ^{n-d^nxdM' (1) 



where n{x) is a three component vector field normalized 
as n • n = 1. The mass scale A is a free parameter and 
in this paper we set A = 1. Stable soliton solutions exist 
when gi > 0. 

The static field n{x) maps n : B? ^ and the con- 
figurations are classified by the topological maps charac- 
terized by a topological invariant H called Hopf charge 



H 
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A^F, F = dA 



(2) 
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where F is the field strength and can be written as F = 
[n ■ dn A dn). 
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The static energy Est t from the action has a topo- 
logical lower bound , 



(3) 



where K = 16n y/gl. 

Performing numerical simulation, one can find that 
the static configurations for H = 1,2 have axial symme- 
try H. Thus "the toroidal ansatz" which was studied in 
Ref.Q is suitable to be imposed on these configurations. 
The ansatz is given by 



ni = y^l — w'^{i],(3) cos{Na + v{r], (i)) , 

n2 = \/l - w2(?7,/3) sin(A^a + vij], (3) , (4) 

where {rj, (3, a) is toroidal coordinates which are related 
to the B? as follows: 



a sinh n cos a 
X = ,y 



a sinh rj sin a 



asin/3 



(5) 



with T = cosh r/ — cos [3. 

The function w{'q, (3) is subject to the boundary con- 
ditions w(0,/3) — f,w(oo,/3) = —1 and is periodic in (3. 
v{ri, P) is set to be ^(77, (3) = M(3 + voir], (3) and t;o(, (3) is 
considered as a constant map. Equation Q then gives 
H = NM. 

In this paper we adopt a simpler ansatz than Q , which 
is defined by 



ni = v/f ~w^{t]) cos{Na + M/3) , 
n2 = Vl - w^iij) sm{Na + Mf3) , 



(6) 



where 10(77) satisfies the boundary conditions w{0) = 
l,w{oo) = —1. We numerically study soliton solutions 
for both ansatz and ((HJ. By comparing those results, 
we found that this simple ansatz produces at most 10% 
errors and it does not affect to the property of the soliton 
solution. 

By using ©, the static energy is written in terms of 
the function w{r]) as 



27r a / dr] 



+ sinh ?7 cosh ry(w')'^[/M,Af (77) 



iw') 



l\2 



w 



— , Um.n{v) = [M^ + --^) ■ 
drj \ smh rj ' 

The Euler-lagrange equation of motion is then derived as 

,/2 



W 



WW 



l-W^ (1 - w2)2 



UM,N{'n)w 



^ ( -2N'^ coth^ rjw' + (cosh^ rj + sinh^ v)UM,Niv)w' 



2a2 

+ sinh r] cosh 'qUM,N {'n)w" ) = . 



(7) 



The variation with respect to a produces the equation for 
variable a. Soliton solutions are obtained by solving the 
equations for a as well as for w. 



III. EFFECTIVE ACTION IN THE 
YANG-MILLS THEORY WITH CFNS 
DECOMPOSITION 

In this section, we briefly review how to derive the 
SFN effective action from the action of SU(2) Yang-Mills 
theory in the infrared limit 12, T^. For the gauge fields 
A^, the Cho-Faddeev-Niemi-Shabanov decomposition is 
applied [3 IH El El 



(8) 



The first two terms are the "electric" and "magnetic" 
Abclian connection, and are chosen so as to orthog- 
onal to n, ■ n — 0. Obviously, the degrees of free- 
dom on the left- and right-hand side of Eq.® do not 
match. While the LHS describes Scoior x 4Lorentz = 12, the 
RHS is comprised of (C^ :)4Lorentz + {n :)2coior + {W^ : 
)3coior X 4Lorentz - 'in-w^.^o = 14 degrees freedom. Sha- 
banov introduced in his paper jl2| the following con- 
straint 



x(n,C^,W^) =0, with x-n = 0. 



(9) 



The generating functional of YM theory can be written 
by using Eq.® as 



Z = J VnVCVWS{x)AFpAse- 



(10) 



App and Sgi are the Faddeev-Popov determinant and 
the gauge fixing action respectively, and Shabanov in- 
troduced another determinant Ag corresponding to the 
condition % = 0. YM and the gauge fixing action is given 
by 



S- 



YM 



gf 



r 1 

Ag 



2 ^l-''^ ' ^1^1^ 



1 



2ag.g' 



Inserting Eq.Q into the action, one obtains the vacuum 
functional 



Z 



(11) 



with 



= -d'^S^iy - dfj,n (g) duTi + dun , 
Q^i. = 9f^ndu + duTbd^ + d^dtjU , (12) 
^^«/ = di,{n ■ diyU X d^n) + d^n ■ d^n x n , 

^ df,{nx d'^n) , (in gauge ag = 1) 



3 



and 



= + Ws ■ QsxM^ll . (13) 



The classical action of n including the gauge fixing term 
is given by 



So 



(14) 



The 5 functional is expressed by its Fourier transform 



(15) 

Integrating over C, W , cf>, wc finally obtain 

g-Sc«- ^ e-^-AFpAs(detAf^)-i/2(detM^)-i/2 
x(det-0^(M^);,iQ^)-i/2 (16) 

where several nonlocal terms and the higher derivative 
components have been neglected. 

We perform the derivative expansion for the four de- 
terminants in Eq. (|16|l under the following assumptions 

(i) the theory is valid for the momenta p with k < p < 
A (fc, A are infrared and ultraviolet cut-off) 

(h) \dn\ < k 

(iii) the higher derivative terms, such as d^n are omit- 
ted. 



The effective action is then given by 
5'cff = / d'^x 



^(W + |(a^nxa,n)2 



(17) 



For gi > and g2 — 0, the action is identical to the FSN 
effective action 

In order to get the stable soliton solutions, g2 must be 
positive 0. However, g2 is found to be negative accord- 
ing to the above analysis. Therefore we consider higher- 
derivative terms and investigate if the model with the 
higher-derivatives possess soliton solutions. 



IV. SEARCH FOR THE STABLE SOLITON 
SOLUTIONS (1) 

The static energy is derived from Ea. (|17|l as 
rl 



d-'x 



E2{n) + Ei'\n) + Ei''\n) 



(18) 



A spatial scaling behaviour of the static energy, so called 
Derrick's scaling argument, can be applied to examine 
the stability of the soliton 19]. Considering the map 



X 1-^ X' 

energy scales as 



fix (/i > 0), with n 



n{fix), the static 



E2 (n(^) ) + 4'' ("*^^ ) + (n(^) ) 



= -E2{n)+^iiEi'\n)+El'>in)). (19) 

Derrick's theorem states that if the function e(/i) has no 
stationary point, the theory has no static solutions of the 
field equation with finite density, other than the vacuum. 
Conversely, if e(/i) has stationary point, the possibility 
of having finite energy soliton solutions is not excluded. 

Eq.lO is stationary at ^ = \/e2/{e[^^ +£'f^)- Then, 
the following inequality 

gi{din X djU^ + g2{dinf{djnf 

= g^{d,nf{d,nf - g^{d,n ■ d,nf + g2{d.,nf{d,nf 



(2)/ 



> 



g2{^^n ■ d,nf (•.■ {d,nf{djnf ^ {^^n ■ d^nf) 



(20) 



ensures the possibility of existence of the stable soliton 
solutions for 92 = 0. As mentioned in the section UTTl g2 
should be negative at least within our derivative expan- 
sion analysis of YM theory. 

A promising idea to tackle the problem was suggested 
by Gies • He considered the following type of effective 
action, accompanying second derivative term 



'5'cff — 



dS 



i(a^n)2-f |-(a^nxa,n)2 



(21) 



Here we choose positive value of g2 and assign the explicit 
negative sign to the third term. In principle, it is possible 
to estimate the second derivative term by the derivative 
expansion without neglecting throughout the calculation. 
The calculation is, however, very laborious and hence we 
show only one simple example of the C determinant. The 
determinant is real and thus it is expanded as follows 

log(detM^)-i/2 = _lTrlog(-92 + Q^n ■ d^n) 



--Trlog(a'*) 



Trlog 



(9n2) (an)4 d^{dnf 



92 



94 



54 



Trlog(9^) 



1. 



-Tr 



Q2 



Q2 



94 

+ 0{{dnf) 



94 



(22) 



4 



where we have defined d^n-d^n (5n)^, {d^n-d^rif 
{dnY . Employing the integral formulas (Ijj 



d^p 1 



1 



[fe,A] (2^)4 p2 



167r2 
1 



(A2 - fc2) 



together with the equality (up to second derivative) 

d^{dt,n ■ 9^n) = -d^{n ■ d^n) = -d^n ■ d^n , (23) 
we obtain the form for the C determinant 



log(detA'/^)"i/2 ^ 
1 



327r2 



A 



(A^ - k'){d^ny - log -(a^n X d,nY 



A, 



A 



log^(9^n)'* - log— (9^n • d'^n) 



(24) 



The other determinants can be estimated in a similar 
manner. 

The static energy of Eq.|2U with the ansatz is 
written as 



Estt = 27r^a / drj 
, 91 



sinh T] cosh 7]{w')'^U m.n (v) 



+ {l-w^)UM,Niv) 



4a2 

- sinh n cosh n 



\ ~ w 



2 +{1-w^)Um,n{v) 



92 

'4a2 



^coth 1] + sinh^ r; — sinh r] cosh rj^ - 
+ (sinh J] cosh 77 — sinh^ 77) ( 1 — w^)M^ 

ITT 



+2 



sinh 77 cosh 77 1 — 



1 — 



+ (l-7«2)^^,^^(^)]' + (^")2| 



d 7« 



The Euler-Lagrange equation of motion is derived by 

= 0, (25) 



d^^dEstt')^ ^ d_/dEstt_\ dEstt 
dri\ dw' J 



drf' V dw" 



dw 



which is too complicated and thus we adopt the following 
notation 

fo{w,w',w") + gifi{w,w' ,w") 

+g2f2{w,w' ,w" ,w^^\w^'^^) = . (26) 

Here w*-^-* , w^'*-' represent the third and the fourth deriva- 
tive with respect to 77. The first two terms of Ea. (|25|l are 
identical to those in Eq.JTjl. 



Unfortunately, we could not find out stable soliton so- 
lutions from Eq. H26|l for any value of 52 ■ 
From the relation 



d^x[{d^n ■ d^n) - (d^n)^] = / d^x{d^n x n)^ ,(27) 



one easily finds that the static energy obtained from the 
last two terms in Ea. (|21|l 



d^xid'^n X nf 



(28) 



gives the positive contribution. The total static energy is 

stationary at ^ \J E2/{e[^^ + E^-') and hence the pos- 
sibility of existence of soliton solutions is not excluded. 
And also, the positivity of Ea. (|28|) does not spoil the 
lower bound Q of original SFN and the possibility is 
still not excluded, too. 

Therefore, we suspect that the absence of the stable 
soliton is caused by the fact that higher derivative theory 
has no lower bound state. We shall investigate the lower 
bound in the higher derivative theory in detail in the next 
section. 



V. HIGHER DERIVATIVE THEORY 

In this section, we address the basic problems in the 
higher derivative theory 0, |2^ [2^ |3 US which essen- 
tially falls into two categories. The first problem concerns 
the increase in the number of degrees of freedom. For 
example, if the theory contains second derivative terms, 
the equation of motion becomes up to the order in the 
fourth derivative. Thus, four parameters are required for 
the initial conditions. If one considers more higher order 
terms, the situation gets worse. However, this is not se- 
rious problem for our study because our concern is the 
existence of static soliton solutions. The second problem 
is that the actions of the theory are not bounded from 
below. This feature makes the higher derivative theories 
unstable. 

The lagrangian and the hamiltonian formalism with 
higher derivative was firstly developed by Ostrograd- 
ski |23|. We consider the lagrangian containing up to 
nth order derivatives 



I dtCiq,q,--- ,q^^^) 



(29) 



Taking the variation of the action 5S — leads the Euler- 
lagrange equation of motion 



(30) 



The hamiltonian is obtained by introducing n generalized 
momenta 



EM) 



j=i+i 



dt^-'-^ \dq 



5 



dC 



Pi = TTTT^ - —Pi+i ^ i=l,---,n-l, (32) 



dt 

and n independent variables 



91 = 9, 

g, = q(*-i), (i = 2,... ,n). 



(33) 



The lagrangian now depends on the n coordinates qi and 
on the first derivative g„ — q^^^K The hamiltonian is 
defined as 

n n— 1 

H{qi,Pi) = ^Piqi - £ = ^piqi+i +p„9„ - C . (34) 

1=1 i=l 

The canonical equations of motion turn out to be 



dn 



dn 



Qi n 1 Pi n 

opi oqi 



(35) 



We consider a simple example including second deriva- 
tive term ||25j, defined as 

/: = i(l + £W-^-V-^eY, (36) 

where constant e works as a coupling constant of second 
derivative term. The equation of motion is 

(1 + eV)q + cj2g + e2^(4) = 0. (37) 

From Eq. 132II . one gets 

dC d /dC\ , 9 9.. 9--- 
Thus the hamiltonian becomes 

H — XTTq + qHq ~ C 

= i^i - - ^(1 + + ^^'9' . (39) 

We introduce the new canonical variables 
1 



9+ 

q- 

p- 



(e w g-TTg) , 



{e UJ q-TTq), 



and the hamiltonian has of the form by using these vari- 
ables 



n^\{pl+u\l)-t(^pl + ±,ql) 



2^' 



The corresponding energy spectra is then given by 

E^{n+]-)LO-{m+]-)- , n,m = 0,1,2,- •• (40) 
2 2 e 

One can see that in the limit e ^ the energy goes to 
negative infinity rather than approaching to the harmonic 
oscillator energy eigenstates. 

To obtain physically meaningful solutions, we employ 
the perturbative analysis where the solution is expanded 
in terms of the small coupling constant and the Euler- 
Lagrange equation of motion is replaced with the cor- 
responding perturbative equation. The solutions of the 
equations of motion that are ill behaved in the limit e — > 
are excluded from the very beginning 23, 2^, 2M ■ 

We assume that the solution of Eq. H37I) can be written 

as 



9pcrt 



(t) = £6>(0. 



(41) 



n=0 



Substituting Eq. (|41|l into Eq. (|37|l and taking time deriva- 
tives of these equations, we obtain the constraints for 
higher derivative terms 

0(e°) 

equation : qo + uP'qa = , (42) 
constraints : Q = —uj'^qo, '4 = w^9o ■ (43) 
O(e^) 

equation : q2 + uj'^qo + ^^^92 + 9 = 0, 

q2+Lu'^q2 = 0, (using gai, 63) ,(44) 
constraints : 4 2 = ^^^^92, 9 2 = ^^'*92 ■ (45) 

equation : q^ + uP'q2 + ijJ^qi -f 9 2 = 0, 

=^ 94 + ^^^94 = , (using (gSl) ,(46) 
constraints : q 4 — —uj'^qi, '4 4^ — UJ^q^ . (47) 

Combining these results, we find the perturbative equa- 
tion of motion up to 0{e*) 

9pert+t^^9port =0(e^). (48) 

which is the equation for harmonic oscillator. 



VI. SEARCH FOR THE STABLE SOLITON 
SOLUTIONS (2) - PERTURBATIVE 
ANALYSIS - 

Let us employ the perturbative method introduced in 
the last section to our problem. We assume that 172 is 
relatively small and can be considered as a perturba- 
tive coupling constant. Thus, the perturbative solution 
is written by a power series in g2 



w{r]) 



(49) 



n=0 



6 
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FIG. 1: The function w{-q) for gi = 0.4, 32 = 0, 0.05, 0.1 (the 
rescahng radial coordinate x = r]/{l ~ ij) is used). 
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FIG. 2: The energy density for gi = 0.4, g2 = 0, 0.05, 0.1. 

Substituting Ea. (|49|l into Ea. (|26|l . we obtain the classical 
field equation in 0(52) 

/o(wo, w'q, Wq ) + 51/1 (wo, w'q, Wq ) = . (50) 

Taking derivatives for both sides in Ea. H5U|) and solving 
for Wg , Wq^' , Wq^' read the following form of the constraint 
equations for higher derivatives 

w^^ = F'^'\wo,w'^) , z = 2,3,4. (51) 
The equation in 0{g\) can be written as 

(/o +.92/1)0(9!) + f2{w(),w'Q,WQ,wf\w\^'^) = 0. (52) 

Substituting the constraint equations (|51|l into Ea. H52|) 
and eliminate the higher derivative terms, one can obtain 



224-1 . 1 . 1 . 1 . 1 . 1 . 1 . 1 . 

222- 

■ 

220- 

218- , 
216- 

214- . . - 

212- . * 

210-, 

208- 

206^^ • ^ • ^ • 1 • 1 • ^ • ^ • 1 • 

0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 



FIG. 3: The energy as a function of (?2((7i = 0.4). 

the perturbative equation of motion 

Uw, w\ w") + g^h{w, w',w") + g2h{w. w') = 0{gl) . 

(53) 

Now Ea. (|53() has topological soliton solutions. Our re- 
sults of the estimated function w{ri) and the energy den- 
sity are displayed in Figs lll2l (In all figures, we show the 
results for the case of Hopf charge iJ = 2;A^ = 2,M = l). 
We have small changes for varying the coupling constant 
52 • The dependence of the (72 for the total energy is shown 
in Fig|3| It can be seen that the change is moderate with 
respect to 172 • 

VII. SUMMARY 

In this paper we have studied the Skyrme-Faddeev- 
Niemi model and its extensions by introducing the reduc- 
tion scheme of the SU(2) Yang-Mills theory to the corre- 
sponding low-energy effective model. The requirement of 
consistency between the low-energy effecive action of the 
YM and the SEN type model lead us to take into account 
second derivative terms in the action. However, we found 
that such an action including the second derivative terms 
does not have stable soliton solutions. This is due to the 
absence of the energy bound in higher derivative theory. 
This fact inspired us to employ the perturbative analysis 
to our effective action. Within the perturbative analysis, 
we were able to obtain the topological soliton solutions. 

Our analysis is based on perturbation and the coupling 
constant 172 is assumed to be small. However, Wilso- 
nian renormalization analysis of YM theory [3| suggest 
that the coupling constants gi , 32 (and the mass scale 
parameter A) depend on the renormalization group time 
t — log/c/A (fc,A are infrared, ultraviolet cutoff param- 
eter) and those are almost comparable. To improve the 



7 



analysis, we could perform the next order of perturba- 
tion, but it is tedius and spoils the simplicity of the FSN 
model unfortunately. 

It should be noted that our solutions do not much dif- 
fer from the solution of original SFN model, at least in 
the perturbative regime. We suspect that some appro- 
priate truncation {like "extra fourth order term + second 
derivative term" ) always supply the stable solutions that 
are close to the original SFN model. Thus we conclude 
that the topological soliton model comprised of the "ki- 



netic term + a special fourth order term" like SFN model 
is a good approximation. 
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